Ordering in a frustrated pyrochlore antiferromagnet proximate to a spin liquid 
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We perform a general study of spin ordering on the pyrochlore lattice with a 3:1 proportionality of 
two spin polarizations. Equivalently, this describes valence bond solid conformations of a quantum 
dimer model on the diamond lattice. We determine the set of likely low temperature ordered phases, 
on the assumption that the ordering is weak, i.e the system is close to a "U(l)" quantum spin liquid 
in which the 3:1 proportionality is maintained but the spins are strongly fluctuating. The nature 
of the 9 ordered states we find is determined by a "projective symmetry" analysis. All the phases 
exhibit translational and rotational symmetry breaking, with an enlarged unit cell containing 4 to 64 
primitive cells of the underlying pyrochlore. The simplest of the 9 phases is the same "R" state found 
earlier in a theoretical study of the ordering on the magnetization plateau in the 5 = 3/2 materials 
CdCr2 04 and HgC^CU- We suggest that the spin/dimer model proposed therein undergoes a direct 
transition from the spin liquid to the R state, and describe a field theory for the universal properties 
of this critical point, at zero and non-zero temperatures. 

PACS numbers: 64.60.-i,71.10.-w,71.27.+a 



I. INTRODUCTION 

Charge and/or magnetic order is an apparently central 
feature of the ground states of Mott insulators. The de- 
tailed nature of this order can be difficult to understand 
theoretically, particularly when "frustration" is present. 
By frustration, we mean the presence of competing in- 
teractions, which lead, in some appropriate "classical" 
limit, to a large number of degenerate ground states. 
This degeneracy is lifted by fluctuations, thermal and 
quantum, or additional interactions beyond those of the 
classical limit. However, the classical degeneracy can be 
lifted in many different ways, making the ultimate ground 
state very sensitive to details of the Hamiltonian. Apart 
from demanding an extremely detailed microscopic un- 
derstanding of a given material (always desirable, but 
not so easy to come by!), is there any way to attack such 
problems? The approach we follow in this paper is to pre- 
sume that fluctuations (in this case quantum) are strong, 
which requires that the ordering itself is weak (i.e. the 
charge and/or spin modulations are small in amplitude). 
If so, we may presume the system to be close to some "liq- 
uid" state, in which no order (in the conventional sense 
- but see below) is present. One may then explore the 
possible ordered states which occur as weak instabilities 
of the liquid. 

Recently, such a view to charge ordered states of two- 
dimensional lattice boson systems has been systemati- 
cally pur sued in Refs. [l||2j, with specific precedents in 
Refs. I-I ll In that case, the liquid state was taken to 
be a superfluid. There, the possible charge ordered Mott 
insulating states proximate to the superfluid were dis- 
cussed by considering the instabilities due to proliferation 
of vortices. This was made systematic by uncovering the 
multiplet structure of the vortex states, determined by 
symmetry. In particular, vortices were shown to trans- 
form under a projective representation of the lattice space 



group, or Projective Space/Symmetry Group (PSG)^&. 
The PSG was shown to determine the structure of the 
action of the critical theory for the superfluid-Mott tran- 
sition, and the nature of the possible charge ordered Mott 
phases. The PSG depends upon the lattice symmetries 
(space group) and the mean conserved boson density. All 
the above considerations for bosons apply equally well to 
spin models with U(l) rather than full SU(2) symmetry, 
with the conserved S z taking the role of boson "charge" . 
This situation is not uncommon, as it is realized when- 
ever an approximately isotropic magnet is subjected to a 
uniform magnetic (Zccman) field. We will focus on this 
realization here. 

In this paper, we will apply an analogous set of ideas 
to bosons/spins on the three-dimensional pyrochlore lat- 
tice. The pyrochlore lattice, consisting of corner-sharing 
tetrahedra (Fig. , takes a central role in the study of 
geometrical frustration in three dimensions. A number of 
materials, in which electronic and/or spin degrees of free- 
dom reside on this lattice, have been intensely studied in 
this light in recent years. Theoretically, the Heisenberg 
antiferromagnet on the pyrochlore lattice is interesting 
as a candidate "spin liquid" (Ref. |5||y), in which the fluc- 
tuations - thermal or quantum - amongst frustration- 
induced degenerate quasi-ground states prevent the oc- 
currence of long-range magnetic order at temperatures 
well below the Curie- Weiss temperature, possibly all the 
way down to T = 0. Quantum spin liquid states can sus- 
tain anomalous spin- 1/2 spinon excitations^, forbidden 
in conventional phases of matter. 

We therefore choose to take as our proximate liquid 
phase not a superfluid (or magnetically ordered phase), 
but instead a particular quantum spin liquid, a so- 
called "U(l)" spin liquid state (see Ref. |8|and references 
therein). A general feature (e.g. on different varieties of 
lattices) of such states is that they exhibit an emergent 
electromagnetism, i.e. they support an "artificial pho- 
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ton", and excitations can carry emergent U(l) electric 
and magnetic gauge charges (see Ref. kill Oil 1(1 and refer- 
ences therein). The spinons carry the elementary quanta 
of the electric charge (of both positive and negative sign 
in two species of spinon). Another, gapped, topological 
excitation, a monopole, carries the dual magnetic gauge 
charge. A transition out of the spin liquid state to a 
state without broken continuous symmetries is generally 
described as a condensation of these monopoles. The na- 
ture of such a transition, and of the proximate spatially- 
ordered states occuring on the other side of the transition 
is determined by the monopole PSG. The ingredients de- 
termining this PSG, as explored on the cubic lattice in 
Refs llllia are the lattice symmetries and the values of 
some conserved "background" U(l) gauge charges, which 
characterize different U(l) liquid states. 

To fix these background charges, we will focus on a 
specific model containing a U(l) spin liquid phase on the 
pyrochlore lattice. In Refll3l this model was argued to 
describe the physics on the magnetization plateaus ob- 
served recently in CdCr 2 04 , HgCr 2 Q4 14 i 15 , which are 
spin-3/2 antiferromagnets with this lattice structure. In 
particular, the model presumes a local constraint (which 
may be understood as the restriction to the classical 
ground state subspace) of three Cr spins fully polar- 
ized (S z = +3/2) along the applied field, and another 
fully polarized antiparallel to it (S z = —3/2), on each 
tetrahedrorM^. We refer to this condition as the "3:1" 
constraint. The model of Ref^^ arises as an effective 
Hamiltonian in this constrained subspace, and takes the 
approximate form: 

Wqdm = V^2(\oa)(o a \ + \o b )(ob\) 
p 

-Kj2(\°A)(o B \ + h.c), (1) 
p 

where indicates a sum over all hexagonal plaque- 
ttes on the pyrochlore lattice, and \Oa), |Ob) are spe- 
cific states with alternating majority and minority spins 
(ITITITI) and IITITIT)) on the given plaquette. This 
model is exactly equivalent to a number of other mod- 
els in the theoretical literature. First, it can be mapped 
directly to a quantum dimer model on the diamond lat- 
tice, the diamond lattice sites being centers of pyrochlore 
tetrahedra (see Sec.[nJ). A number of such dimer models 
have been considered in the literaturei£ii£*i2*22i2L22. Sec- 
ond, the dimer model in turn can be rewritten as a par- 
ticular compact U(l) gauge theory. The 3:1 constraint of 
the spin model maps directly to the background charge 
of this gauge theory. In this way, the essential ingredients 
fixing the monopole PSG are determined. A systematic 
analysis of the spatially-ordered states proximate to the 
spin liquid is therefore possible, and is the main subject 
pursued in this paper. 

More microscopically, it is possible to show that the 
spin/dimer model of Eq. Q indeed exhibits a U(l) spin 
liquid ground state when the dimensionless parameter 




FIG. 1: (Color online.) The pyrochlore lattice structures, 
shown as a network of corner sharing tetrahedra. The atoms 
occupy the corners of the teterahedra. 



v = V/K satisfies v c < v < 1. This argument, analo- 
gous to the ones in RefsHH is described in Sec. IIIBI 
The critical coupling v c is not known, but based on nu- 
merical analysis of other similar models probably satis- 
fies v c > —0.5 or so 20 i 21 i 22 i 23 . For the application to 
CdCr204 , HgCr204 , it was estimated in Ref 1 131 that 
v ps —1.2 < v c . The nature of the ground state in that 
case may perhaps be more accurately understood by ex- 
trapolation from the limit v — * — oo. The ground state 
can be determined classically in that limit, and in Refll3l 
was found to exhibit a particular spatial ordering pattern 
with a quadrupled unit cell. It can be understood by 
order-by-disordeB 2 ^ 2 ^^ reasoning as the classical state 
with the most possible "resonances" - off-diagonal quan- 
tum moves via the K term in Eq. Q - to other state. 
We therefore refer to it as the "R" state. 

The alternate approach, which we pursue here, is to 
approach the physical limit from the spin liquid state, 
asking which ordering pattern emerges from the PSG 
analysis. Remarkably, we find that the simplest possi- 
ble ordered state proximate to the U(l ) spin liquid is the 
R state! This suggests the possibility that the R state 
in the physical limit may be close to a phase transition 
to the spin liquid state. The analysis of this paper pro- 
vides an analytical framework for such a transition, both 
of quantum and of thermal nature. 

Apart from the simplest R state, a number of or- 
dered phases come out of this analysis, and are shown 
in Figs. 1111191 A salient feature of all these phases is 
that they exhibit an enlarged unit cell relative to that of 
the original pyrochlore lattice. Interestingly, this set of 
phases therefore does not contain the simplest "ferrimag- 
netic" ordered state illustrated in Fig. [3 in which the unit 
cell is not enlarged though the point group symmetry of 
the crystal is lowered. This indicates that the ground 
states selected out of the classically degenerate ground 
state manifold are identifiably different from those pre- 
ferred by other degeneracy breaking mechanisms, such 
as e.g antiferromagnetic second-neighbor exchange cou- 
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FIG. 2: The simplest ordered state consistent with the 3:1 
proportion of majority to minority sites, which does not ex- 
hibit an enlarged unit cell, and has only a four-fold degener- 
acy. In the spin language, this state has three spins aligned 
with the field and one anti-aligned, so we will denote this as 
a "ferrimagnetic" state, analogous to the ferrimagnetic state 
found in two-dimensional triangular lattice antiferromagnets 
in a field. This state is not generically proximate to the U(l) 
spin liquid. 



pling. Identification of the precise nature of the ordered 
state in experiments therefore indirectly gives useful 
formation on the importance of quantum fluctuations 



m- 



The remainder of the paper is organized in the follow- 
ing way. In Sec. [H] we describe our theoretical model 
(lattice QED). In Sec. HHI we illustrate and motivate the 
dual transformation we make on the Hamiltonian intro- 
duced in the previous section. In Sec. II VI we derive and 
present the effective action for the monopole degrees of 
freedom that appear in the dual theory. We then study 
the action in the mean-field approximation in Sec.lVland 
present the resulting charge ordered phases in Sec. I VII In 
Sec. IVIII we carry out a renormalization group analysis 



of the action. Finally in Sec. IVIIII we conclude with a 
discussion of the phase diagram and critical phenomena 
in the spin/dimer model at non-zero temperatures. Im- 
portant but lengthy formulas and results are given in the 
appendices. 



A. Equivalence to compact lattice QED 

First, it is useful to discuss how the model can be cast 
into a lattice U(l) gauge theory. A reader not famil- 
iar with lattice QED may wish to consult the review by 
Kogut, Ref. |23- As mentioned in the introduction, it is 
first convenient to pass from the pyrochlore to the dia- 
mond lattice. This is accomplished by focusing on the 
centers of the tetrahedra (labeled by the sites a and b in 
Fig. |3| that make up the pyrochlore. The centers of these 
tetrahedra make up a diamond lattice. Each site on the 
pyrochlore lattice connects two nearest neighbor tetrahe- 
dra, and can be identified with a link between the centers 
of the two tetrahedra. Thinking in terms of the centers 
of the tetrahedra as the sites of a new lattice, the spins 
sit on the links of a diamond lattice. The spin states may 
therefore be regarded as dimer coverings of the links of 
the diamond lattice, and the effective Hamiltonian as a 
quantum dimer model. 




FIG. 3: Section of a pyrochlore lattice. The pyrochlore sites 
are denoted by i and the tetrahedra are identified by a, b la- 
bels. Drawing links between the tetrahedron centers forms 
the links of a diamond lattice, with the tetrahedron centers 
corresponding to the sites of the diamond lattice. The figure 
indicates the bipartite nature of the diamond lattice, which 
is evident in the notion of up/down pointing tetrahedra as 
indicated by the labels u/d. 



II. THEORETICAL MODEL: FROM SPINS TO 
COMPACT LATTICE QED 

In this section, we reformulate the 3:1 quantum 
spin/dimer model of Eq. (JIJ as a lattice U(l) gauge the- 
ory, and describe its phase diagram, which contains both 
spin liquid and ordered states. This Hamiltonian has a 
single dimensionless parameter, v = V/K, so the zero 
temperature phase diagram is entirely determined by v 
(we fix K > by convention - its sign can be changed by 
a suitable canonical transformation, and has no signifi- 
cance) . We discuss the structure of this phase diagram. 



The gauge nature of the problem is simply a conse- 
quence of the local 3:1 constraint on each pyrochlore 
tetrahedron, or equivalently, that each diamond site is 
covered by a single dimer. To map this onto a more con- 
ventional gauge theory, we note further that the diamond 
lattice is bipartite, so we can define "up" (it) and "down" 
(d) sub-lattices corresponding to neighboring tetrahedra 
in the original pyrochlore lattice, as illustrated in Fig. [3] 
We can define thereby a discrete oriented electric field 
variable on the diamond lattice links, equal to zero if the 
corresponding pyrochlore site is a majority site, and equal 
to ±1 if the pyrochlore site is a minority one, choosing the 
field to always point from the up to the down diamond 
sublattice. We specify the spin/dimer configurations on 
the pyrochlore by a discrete variable hi (i is a pyrochlore 
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lattice site), such that 



majority site, 

1 minority site. 



(2) 



Mathematically, for the site i lying on the diamond 
link ab, 



where 



E a b — £a,n a b j 



+1 , a 6 u 
— 1 , a e d 



(3) 



(4) 



The electric field direction can be identified by the index 
ordering that gives it a positive value, as in Fig. 

The model we construct requires that each tetrahedron 
must include 1 minority site, 



N„ 



l 



(5) 



where the label a identifies the various tetrahedra, and % 
is summed over pyrochlore lattice sites on tetrahedron a. 
See Fig. [31 for an illustration. This 3:1 constraint on the 
pyrochlore lattice in the QED formulation maps directly 
into 



divE , 



(6) 



which is rather reminiscent of Gauss' law, where we have 
used the lattice divergence div E — J2b Eab- The "charge 
distribution" in this picture is that of alternating positive 
and negative charges on the diamond lattice. Positive 
charges sit on the u sub-lattice, and negative charges sit 
on the d sub-lattice. 

A further consequence of this mapping is the presence 
of global topological charges, which are conserved in pe- 
riodic (or infinite) systems. In particular, if one draws 
any surface not passing directly through diamond sites, 
the net electric flux through this surface is conserved. 
If this surface is compact and closed, this gives no ad- 
ditional information beyond the Gauss' law of Eq.©. 
However, one may also consider non-compact surfaces 
that extend across the entire sample, and the electric 
flux through such surfaces is not determined by Gauss' 
law. A sufficient set of surfaces are the four non-parallel 
flat planes containing two-dimensional triangular lattices 
of pyrochlore sites (any two such parallel planes have 
the same electric flux). We denote the corresponding 
four electric fluxes by a four- vector {E±, E2, E3, E4). The 
fluxes may be chosen positive, and can in principle take 
any integer value from < Ei < where Aa is 

the number of triangular sites in the plane. By the 3:1 
(Gauss' law) constraint, the sum E\ + E 2 + E 3 + E4 = Aa 
is fixed. The electric flux sector containing the ground 
state varies with v. 

It is conceptually useful to also map the full Hamil- 
tonian of the model to a form more familiar in lattice 



QED. To do so, we must introduce the phase operator . 
conjugate to the number operator hi, satisfying 



4>j, fii 



-i 8j 



(7) 



where Sij is the Kronecker delta function. The operator 
e + l 4>j creates a minority site at site j, and using it we can 
construct any hopping term we wish for minority sites 
(down spins). With the canonical "rotor" variables hi, 
and 4>i, in principle an infinite set of number states with 
all integer eigenvalues of hi are allowed. To faithfully 
represent the original spin/dimer model, therefore, we 
will include a large term U in the Hamiltonian which, 
in the limit U — > 00, restricts the site occupancies to 
hi = 0,1 as desired. We thereby obtain 

H =| Yl * - J ) + U * - l ? 

i a 

+ V *^2 I/" 1 ' i'WckW 1^4,0^5,1^6,0 + (n» <-> 1 - rii) 



o 



K 
~2 



o 



(8) 



Here O denotes the hexagonal plaquettes on the py- 
rochlore lattice, and the indices on <f> enumerate the site 
(links) on the hexagon. The constraint operators A Q 
commute with TL by construction, so for sufficiently large 
Ut the ground state will indeed satisfy Eq JSJ. Moreover, 
when the constraint is enforced, the Ut term plays no fur- 
ther role. Formally, we are principally interested in the 
limit U/K — > 00, as described above. 

The K term when rewritten in this way appears as 
a rather complicated-looking multi-particle hopping am- 
plitude. In fact, this form is actually the simplest one 
allowed by the constraint (|3J). The hopping of a down 
spin from one lattice site to another can be decomposed 
into a series of hops along nearest neighbor links, so it is 
sufficient to analyze the simplest allowed moves. In gen- 
eral, the hopping on a nearest-neighbor link will violate 
O on 2 separate tetrahedra. Thus, any series of such 
hopping events will do the same. Analysis shows that it 
is only possible to hop between tetrahedra along closed 
contours. On the pyrochlore lattice the smallest closed 
contours are hexagonal plaquettes. Any other closed con- 
tour on the lattice can be constructed from these min- 
imal moves, so we shall consider exclusively this "ring 
exchange hopping" on the hexagonal plaquettes. Such 
moves are illustrated in Fig. 0] 

This Hamiltonian can now be re-expressed as a lattice 
gauge theory. Analogously to the definition of the electric 
field, we define the vector potential as 



•Aab 



(9) 



The electromagnetic variables introduced above obey the 
same canonical commutation relations as 6 and h, since 
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(10) 



on the same link, and the commutator is on different 
links. Note that from their definitions, E a b is integer- 
valued (and in particular = 0, ±1) and A a b is a 2n- 

periodic phase variable since the operator e 1 ^ creates a 
particle. 

In the new variables the ring exchange term becomes 



e +i{A 1 +A 2 +A 3 +A4+A 5 +Ae) + c = 

e +ifA-di + h c = 2 cos ( curl _4) Q ( 
where we have introduced the lattice curl: 



(curl .A)o = ^ A?f> 

rr'eO 



(11) 



(12) 



In this form, the previously complicated-looking form of 
the K term becomes transparent. 

After substituting the new variables, the Hamiltonian 
(JSJ) takes the form 



n 



const. 

o 



U 
~2 



(a,b) 



irl£,3 



rl£,-3 



K cos ( cur] .4 

o 



(13) 



where the constant is a result of reorganizing the first 
term of © into a quadratic form. By identifying the 
curl of the vector potential with a magnetic field running 
through the plaquettes, B = curl A, 



n 



const. 



(a.b) 



o 



o 



curl.E, 3 



curlB, 



(14) 



After these manipulations, the spin/dimcr Hamiltonian 
has been formulated as a compact Quantum Electrody- 
namics (QED) in 3+1 dimensions. Eq. I|14|) is very simi- 
lar to the standard form of compact QED, but does dif- 
fer from it by the presence of the rather ugly V term 





FIG. 4: (Color online.) Ring exchange hopping on a hexago- 
nal plaquette. The hopping flips between A and B type pla- 
quettes, and is the only simple hopping event that preserve 
the 1 boson per tetrahedron constraint. 



and the modification of the E 2 term by the e a /2 "back- 
ground field" . Please recall also that we are expected to 
take the limit U — > 00 to recover the spin/dimer model. 
Despite the differences, Eq. 1)140 does share all the same 
internal symmetries as the more conventional QED form. 
It is therefore expected to share the same properties in 
regimes where universality is mandated. 



B. Phase diagram of the quantum spin/dimer 
model 



Let us now return to the microscopic form of the 
spin/dimcr model, and the question of the phase dia- 
gram. We will employ the QED formulation where use- 
ful in this analysis. For v — — > — 00, the off-diagonal 
K term can be neglected, and ground state is deter- 
mined by minimizing the (negative) V term over classical 
spin/dimer configurations. The solution is the R state, 
shown in Fig. ^2 This state has only a discrete degen- 
eracy, and is separated from other excited states by a 
gap (of 0(V)), so it is expected to be stable to pertur- 
bation theory in K. Therefore the R state is the ground 
state for v < v c i, with some v c \ < 1. In the R state, 
the electric flux is equally divided, Ei/N& = 1/4 apart 
from 0(1/N&) corrections for some frustrated boundary 
conditions. 

For v > 1, the ground state can be found by rewriting 
the Hamiltonian as follows: 

Hqdm = l*X(|o A > - |ob»«Oa| - (ObI) 
p 

+ (v - 1)KJ2(\°a)(Oa\ + |Ob)(Ob|) .(15) 
p 

In Eq. 1)15)1 Hqdm has been expressed as a sum of posi- 
tive semi-definite projection operators, with coefficients 
that are all positive for v > 1. Therefore the energy is 
bounded below by zero, and any zero energy state is a 
ground state. In particular, any classical state that con- 
tains no A or B hexagons is automatically a ground state. 
As a simple example, consider the "ferrimagnetic" state, 
with no enlargement of the unit cell but broken rotational 
symmetry. It is described as follows. On the diamond 
lattice, each vertex has 4 links emanating from it, which 
we label by u — 0, 1, 2, 3. The ferrimagnetic state (ori- 
ented along a specific //-direction in space) has the same 
3:1 arrangement on all identical tetrahedra. Explicitly, it 
can be written 



|Ferri) M = \] 



= 1) 



(16) 



where the product is taken over the u sublattice of tetra- 
hedra shown in Fig. EI This state clearly obeys the 3:1 
constraint JSJ - one minority site on the u corner of each 
tetrahedron. And moreover, it contains no A or B plaque- 
ttes, and hence is a zero energy ground state for v > 1. 
The ferrimagnetic state is in the "furthest" topological 



6 



sector from the R state, with E{ = (N&, 0,0,0) (and 
permutations). 

The ground states in this regime are, however, highly 
degenerate, and many other classical configurations are 
possible. All such states are "frozen" in that, for any 
value of v, they are exact eigenstates with trivial dy- 
namics, being annihilated by the off-diagonal K term in 
Eq. p. 

At the point v = 1, Eq. I|15|) simplifies (the last term 
dropping out), but remains the sum of positive semidef- 
inite projectors. This is the so-called Rokhsar-Kivelson 
(RK) point. The ground state space is enlarged, and con- 
tains now states in addition to the frozen ones. In partic- 
ular, many configurations containing A and B hexagons 
are now allowed, provided the first projector in Eq. I|15() 
annihilates the quantum state. One construction of this 
type is especially simple. Take a uniform superposition 
of all possible spin configurations, ignoring the 3:1 con- 
straint. Now project this onto the 3:1 manifold. This is 
known as the RK wavefunction. It can be further bro- 
ken into substates, e.g. by projecting out all the frozen 
states. It can also be projected into any of the electric 
flux sectors. 

On reducing v to values slightly below unity, the second 
term in Eq. (|15l) becomes again non-zero but negative 
semi-definite. This indicates configurations with A and 
B hexagons are now preferred in the ground state. The 
frozen states are then highly excited and energetically 
unfavorable. Instead, by an application of the reasoning 
of Ref0, the ground state can be argued to be a U(l) 
spin liquid state. 

We summarize this argument, and the nature of the 
U(l) state. It is by now well-known that RK wavefunc- 
tions for dimer models on bipartite lattices (like the dia- 
mond) display power-law equal-time correlations. These 
correlations may be understood as arising from the 3:1 
constraint. In particular, the calculation of equal-time 
correlations in the RK state reduces to a problem of three 
dimensional classical statistical mechanics: performing a 
statistical average over discrete electric field configura- 
tions subject to the Gauss' law constraint of Eq. ©. It 
turns out that the long-distance behavior of these corre- 
lations is captured simply by taking an effective classi- 
cal free energy density proportional to \E\ 2 and treating 
£ as a continuous Gaussian (but constrained) variable. 
The resulting correlations have a "dipolar" power-law 
form§i2&. 

By an argument originally due to Henley^, these 
power laws can be understood in the quantum theory as 
follows. Evidently, the discreteness of the electric field is 
unimportant at the RK point. Consequently, it is natural 
to treat E and B as continuous, write an effective action 
quadratic in these fields, and apply the Gauss' law con- 
straint. However, there is an additional feature dictated 
by another peculiarity of the model. As we have pointed 
out already, all possible electric flux sectors are degener- 
ate at the RK point. This implies that there should be 
no cost in the energy to shift E by a uniform constant (at 



the RK point). Henley's argument therefore indicates an 
appropriate effective Hamiltonian density is 

H cS = a\ curl E\ 2 + b\B\ 2 + a(l-v)\E\ 2 . (17) 

The last \E\ 2 term must vanish at the RK point, but is 
expected to become non-zero if one perturbs away from 
it. For v — 1, this form can be shown to precisely repro- 
duce the microscopically calculated correlations of the 
RK wavefunction, with specific constants a, b. For v 
larger than one, the negative \E\ 2 term favors the sec- 
tors with "large electric flux", i.e. the frozen states, as 
expected, and E itself develops a non-zero expectation 
value. For v slightly less than one, the positive \E\ 2 
term instead favors the "minimal electric flux" sector 
with Ei = N A /4. 

At low energies, for v < 1, therefore, it is expected 
that the a term above can be dropped in favor of the 
a term, making the effective Hamiltonian simply that of 
the usual non-compact QED. This indicates the system 
is in the "Coulomb phase" of the gauge theory, which has 
the usual properties expected of 3+1 dimensional electro- 
dynamics. In particular, unit test gauge charges can be 
introduced and interact via bounded 1/r Coulomb po- 
tentials. Such a gauge charge corresponds in the original 
pyrochlore magnet to a "spinon" excitation with frac- 
tional spin ±3/2. Thus the Coulomb phase is indeed a 
spin liquid, and in respect of the U(l) gauge structure, 
this state is called a U(l) spin liquid. 

A key difference from standard "non-compact" QED 
does appear at non-zero energies, as a consequence of 
the fact that the B is defined modulo 2ir. The diver- 
gence of this magnetic field is the sum of magnetic field 
values coming out of the plaquettes enclosing one cell in 
the diamond lattice. With the modulo 2tt redundancy, 
configurations with div B an integer multiple of 2n are 
allowed. Thus, compact QED allows magnetic monopoles 
with quantized magnetic charge. These have a finite en- 
ergy cost and are gapped excitations in the U(l) spin 
liquid. 

Although the Coulomb phase emerges in a non-trivial 
way from the spin/dimer model in the vicinity of the RK 
point, we can mimic its low energy physics more simply. 
In particular, the same phase is obtained from Eq. (|14|> 
by dropping the diagonal V term, and taking non-zero 
but finite U (instead of U — * oo). This gives a "soft- 
ened" model with the same universal properties as the 
original spin/dimer model. It is a remarkable fact that, 
by the preceding arguments, these two "sins" compensate 
themselves and give the proper behavior of the original 
spin/dimer model for v < 1, in the Coulomb phase. 

III. DUALITY AND MONOPOLE 
FORMULATION 

On reducing v from values just below one, eventually 
the spin/dimer model must undergo a transition out of 
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the spin liquid state. The resulting state can be mimicked 
in the softened model by increasing U/K. The eventual 
outcome can be understood as follows. For U ^> K the 
"magnetic field" term is subdominant, and so the electric 
field is a good quantum number. This limit is somewhat 
complicated in Eq. I|14[l. because the U term selects two 
degenerate values E a b = 0, e a as low energy states of each 
bond. Indeed this recovers the original effective spin-3/2 
model, which is of course still non-trivial. However, one 
can readily understand in this limit the basic nature of 
the other phases of the theory. To do so, we imagine gen- 
eralizing the U term to include electric-field interactions 
on nearby bonds. This will generally break the large U 
degeneracy in favor of some particular global arrange- 
ment of E a b — 0, e a values. Because of the discreteness 
of E a b-, deviations from this ground state are likewise dis- 
crete, any local rearrangement of the pattern results in a 
non-zero increase in energy, i.e. there is an energy gap. 
Gauge-neutral excitations are created in this way. Ei- 
ther a set of E a \, fields are modified along links forming a 
closed curve, or a pair of diamond lattice sites (on which 
the Gauss' law constraint is violated) is created, with a 
modified path of E a b fields connecting them. The latter 
corresponds to a particle-antiparticle pair, and costs an 
energy proportional to the length of the path. Hence the 
pair itself is bound, the bound state being gauge neutral. 
The individual gauge-charged excitations are said to be 
confined by the linear potential between them. The con- 
fined phase corresponds to having the n, operators cer- 
tain in (jSJ - a Mott insulating phase, with some sort of 
diagonal ordering. In contrast to the Coulomb phase, the 
B field in the confining phase is very strongly fluctuating, 
so monopoles are no longer good excitations. In fact, it is 
appropriate to think of the confined phase as a Coulomb 
phase which has been destroyed by Bose condensation of 
monopoles. The presence of a delocalized monopole con- 
densate can be thought of as leading to strong (gauge) 
magnetic fluctuations in the ground state. 

We therefore wish to reformulate the lattice QED 
model so that the monopole excitations of the Coulomb 
phase are explicit. Here we follow Hermele et al& with 
slight differences. As the reader will recall, the electric 
and magnetic fields in Maxwell's equations are dual when 
there are no charges or currents present. While there are 
no currents in our system, Eq. JSJ shows we do have a 
charge distribution. The duality transformation is 



ib i 



B = curl .4 



(18) 
(19) 



We have thus introduced an explicit operator for the 
magnetic field B and an "electric vector potential" 
a, whose exponential creates a 'magnetic field' since 
an exponential of A creates an "electric field", via 
EH- Here e( 0) is a classical electric field cre- 
ated by the charge distribution ©, and 



It is convenient to choose to be integer- valued, so 
that curia may also be taken integer-valued. A simple 
choice is to take the classical configuration corresponding 
to one of the 3:1 states, e.g. just e^\+u = e a<^o- 



*-. i 




FIG. 5: (Color online.) Direct and dual diamond lattices are 
dual to one another. The plaquettes of one lattice correspond 
to the links of its dual lattice. 



The Hamiltonian in the dual language takes the form 
U 



H = 



u 

o 



curl< 



,(0) 



K^cosB, (21) 



where we denote the sites of the dual lattice by r, and its 
links by r, r' . The hexagons now denote the plaquettes 
of the dual lattice. (Note the change in the summation 
subscripts on both the first and second terms.) 

The dual fields obey the canonical commutation rela- 
tions 



B r 



0t r 



(22) 



dive^ = divE = e a 



(20) 



and the commutator vanishes for different links. The new 
fields are once again conjugate variables. B r r i is defined 
modulo 2tt - an angular variable, and the a r y variable 
is integer valued. 

Standard manipulations can now be used to "soften" 
the inconvenient integer constraint on a r , r 'i remove the 
periodicity of B r r > , and make the monopole variables ex- 
plicit. The reader is referred to Refs lTfll29l and references 
therein for details. These manipulations are inexact, but 
do not change the structure of the phase diagram in the 
vicinity of the transition from the Coulomb to confining 
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phase. One obtains 



W = ^5>urla-e) 2 + ^£s 2 



o 



(23) 



where now a r , r ' and B are real variables. In Eq. Ij23|) . one 
may freely shift e by a gradient, changing only the overall 
zero of energy, since such a gradient does not couple to 
curl a. We have used this freedom to modify the original 
e (0) + e a /2 terms to 



-8. 



//0 



(24) 



which has no divergence, but has the same curl as the 
original "source" fields. As promised, explicit monopolc 
degrees of freedom have been introduced. A monopolc 
number operator N r is slaved (by a dual Gauss' law con- 
straint) to the B field, 

div.B = 27riV (25) 
It is conjugate to the dual phase 9 r , such that 



N r > 



I 1/ y ' 



(26) 



Eq. (|25|l is another U(l) gauge constraint, so it is not 
surprising that the monopole hopping term w respects a 
dual (non-compact) gauge symmetry. 

Although the monopole number N r (which can be both 
positive or negative reflecting the two signs of flux em- 
anating from a monopole) appears nowhere explicitly in 
Eq. i)23|) . it is implicit through the constraint of Eq. <|25|) . 
In the Coulomb phase for large K, therefore, monopolcs 
are energetically costly (though their energy is finite, as 
is easily verified by integrating the associated B 2 energy 
density), with a gap of O(K). Through the w term, how- 
ever, monopoles do not reside in localized states with 
N r = ±1, but instead in superpositions of such states, 
with only ^ r N r — ±1. As K is decreased, the monopole 
energy gap decreases, and at some point it will reach zero. 
This point corresponds to the confinement transition dis- 
cussed in the previous section. 



IV. FORMALISM - MONOPOLE DEFECTS ON 
THE DIAMOND LATTICE 

A. Ground state manifold 

To understand the confinement transition, we must 
understand the nature of the lowest energy monopole 
and anti- monopole states, which condense at the tran- 
sition. They are equivalent by B — > — B symmetry, so 
it is sufficient to study just the monopole states. The 



ultimate field theory will consist of a relativistic field for 
each member of the monopole multiplet, since the rel- 
ativistic description includes particles and antiparticles 
(here antimonopolcs) on equal footing. We will apply 
for the most part a (dual) mean-field approach, taking 
curl a = e in IL'-it , and neglecting fluctuations of a around 
this value. This is sufficient to analyze the spectrum of 
ordered phases near the U(l) quantum liquid. Fluctua- 
tions will be restored later in Sec. IVIII 

We may thus consider the manifold of states with one 
monopole, i.e. N r = 1 on one and only one site of the 
dual lattice, and N r = on all other sites. Through 
the w term in l|23|L the wavefunction of the monopole 
delocalizes, and is described by a tight-binding model, 
which we may write as 



Hti = -w^2{^(r f )il>(r)( 



h.c.} 



(27) 



(r,r'> 



where (r, r') denotes a summation over nearest neigh- 
bors on the diamond lattice. Here ipj.,ip r are cre- 
ation/annihilation operators for the monopole. Note we 
have absorbed a factor of 2-7T into the vector potential 
relative to (|23|l to make our notation more conventional. 
By our mean-field assumption, a r ^ is a c-number vector 
potential carrying a "flux" (actually electric flux) given 
by life. Since it appears only in a periodic exponential, 
the form in Eq. (|24(1 is equivalent to a flux of 27r/4 = 7r/2 
through each dual plaquette. 

This can be understood as follows. The original com- 
pact QED theory had a staggered background charge of 
e a on each direct lattice site. The monopoles, as magnetic 
charges, see this in the same way electric charges would 
see a staggered lattice of magnetic monopoles and anti- 
monopoles. These 'monopoles' and 'anti-monopoles' (we 
use single quotes to denote the dual view, since these are 
actually the background gauge charges) are distributed 
in an alternating fashion at the center of each cell of the 
dual diamond lattice. The neighboring cells to a cell con- 
taining a 'monopole' all contain 'anti-monopoles'. Each 
'monopole' has a "charge" of 2tt. Since all lattice di- 
rections are equivalent, the 'magnetic' flux going out of 
each face of the cell must be the same, as illustrated in 
Fig. The structure of the diamond lattice is made of 
cells where each cell has 4 faces, as opposed to the cu- 
bic lattice which has 6 faces for each of its cells (cubes). 
Thus we conclude that each face in the diamond lattice 
has a flux of =?■ = \ going through it in the direction 
from a 'monopole' cell to an 'anti-monopole' cell. The 
ip r monopole particle thus experiences Aharonov-Bohm 
fluxes of precisely this sort as it moves through the lat- 
tice. 

It proves convenient to describe the links of the dia- 
mond lattice by (r,r') = (a, ji) where a £ u denote the 
sites of the u sublattice, and fi e {0,1,2,3} enumerate 
the 4 links emanating from each u site. Furthermore, we 
can enumerate the u sublattice sites by r = Y^j=i n j^j 
where Sj [Si = %{y + z), a 2 = f (x + z), a 3 = | (x + y)] 
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FIG. 6: (Color online.) Alternating charge distribution emits 
'magnetic' field lines through the faces of a diamond lattice 
cell. 

are the primitive Bravais lattice vectors of the FCC lat- 
tice, and rij span the integer numbers. We refer to this 
coordinate system as "index" space. 

We shall now focus our attention on finding the ground 
state manifold of this Hamiltonian. First we must hnd an 
appropriate choice of the vector potential giving the de- 
sired flux pattern through the faces inside the lattice. To 
this end, the "index space" notation proves particularly 
useful. One such possible vector potential is 

a (n) = , 

( ~ \ ( 28 ) 
a(n) = (a 1 (ri),a 2 (n),a 3 (ri)) = elQ ■ n\ , 

where Q = § (1,0,-1) and e= (1,1,2). 

We proceed to diagonalize the hopping term. General 
eigenstates cannot be found analytically, however min- 
imum energy eigenstates can. We find 8 ground state 
eigenmodes, denoted <& v and $„ where the indices run 
through v = 0,1,2,3. The details of the 8 eigenmodes 
are are left for Appendix lAl 



the Hamiltonian with the same energy. Such a structure 
therefore constrains the minimum size of ground state 
multiplets. This idea should be familiar from the repre- 
sentation theory of SU(2) in standard quantum mechan- 
ics. 

While the symmetry group of the diamond lattice is the 
space group Fd3m2£, the dual Hamiltonian itself trans- 
forms under a projective symmetry group (PSG), since 
a vector potential a appears explicitly. The PSG dif- 
fers from the space symmetry group by a specific gauge 
transformation accompanying every space symmetry op- 
eration. 

We can describe the entire symmetry group using a 
reduced set of operators - generators - from which any 
symmetry operation can be constructed as a product of 
members in the reduced set. One can find a minimal 
such reduced set. We shall consider the minimal set of 
the point symmetry group, as well as the 3 primitive 
translations. We leave the complete description of the 
PSG to AppendixlBl and only state here that the reduced 
set comprises 

• 3 (FCC) translations tj, 

• a reflection symmetry i, 

• an inversion symmetry p, 

• a 2-fold rotation symmetry ri, 

• a 3-fold rotation symmetry r&. 

Our convention is to denote symmetry operations in the 
space symmetry group by lowercase letters, and the cor- 
responding PSG operations by uppercase letters. The 
space group obeys some group algebra, which is just a 
multiplication table. It can be constructed by imposing 




B. Symmetries 

The 8 k-space minima can be related by the symme- 
tries of the Hamiltonian. Each one of the symmetries 
is represented by an operator that commutes with the 
Hamiltonian. However, some of the symmetry operators 
do not commute with each other. If we choose to repre- 
sent the Hilbert subspace of the minimum energy states 
of the Hamiltonian using a basis of eigenstates that are 
also eigenstates of a symmetry operator Si, then for a 
different operator S2 with which it does not commute 
[SijSy 7^ the basis states will not be eigenstates of 
§2- Thus by acting with S2 on the eigenstates of Si, one 
may generate another linearly independent eigenstate of 
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FIG. 7: (Color online 
tor potential pattern 



Projected diamond plane view of vec- 
links between the honeycomb planes 
have ao — 0. All links with the same vector potential value 
are in the same color. 
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a set of algebraic rules on the generators alone, of the 
form 

li...s„ = l, (29) 

where Sj are generators. The details of these relations 
are left to Appendix O 

To obtain the PSG we must add gauge transforma- 
tions to each of these symmetry operations. Consider a 
symmetry operation S with the following action on the 
lattice sites in real space coordinates s : r — > r'. Now 
add a gauge transformation to accompany the lattice site 
transformation 

S : ip{r) — ► if;(r')e- lA ^ . (30) 

Let us examine what this transformation does to a 
generic hopping term in the Hamiltonian 

S: V t (r 2 )^(ri)e- M - 1 ^ 

(31) 

^{r' 2 )^{r' 1 )e +i { K ^- K ^- a ^) . 

We require that the Hamiltonian be invariant under this 
transformation, and so we must have 

e +i(A(r a )-A(n)-a Pll r a ) = g-^V^ ; (33) 

or put more simply 

A(r 2 ) - A(r x ) - a TuT2 = -a r ^ r ' 2 (mod 2ir) . (33) 

Using this procedure we can find the appropriate gauge 
transformations for each of the symmetry operations in 
our reduced set. We leave the details to Appendix FBI 

The PSG has a modified group algebra. The relations 
among the generators differ only slightly from those in 
the space group, which are given in i|Cl(l . Consider any 
succession of symmetry operations that takes every site 
on the lattice back to itself. The same succession of the 
"gauged" operations can therefore only perform a gauge 
transformation on the lattice model. Since all the gauged 
symmetry operations leave the Hamiltonian invariant, 
and the various terms in the Hamiltonian depend explic- 
itly on the vector potential, any gauge transformation 
that is not uniform (a global gauge transformation) will 
modify the vector potential by introducing non-zero gra- 
dients in the gauge transformation. We conclude, there- 
fore, that any succession of symmetry operations that 
leave the lattice sites in place can only undergo an addi- 
tional global gauge transformation when those operations 
are gauged. A rule in the algebra (|2*5|l is in general mod- 
ified into 

Si...S n = e +lB , (34) 

where 9 is some angle. The PSG algebra differs from the 
space symmetry group algebra (|C1|) only in the following 
rules 

h+i ■ h ■ • *T X = 1 > ( 35 ) 



in the space algebra, and 

T J+1 ■ % ■ fr + \ . ff 1 = -i , (36) 

in the PSG algebra. 

The 8 ground state modes we found span the ground 
state manifold, and so any linear combination of these 
states is also a ground state, 

3 

* = ^(^«*,+^,d^), (37) 

v=0 

where 4>v,a now denote the complex amplitudes of the 8 
eigenmodes. Given the symmetry operations in momen- 
tum space, we find the ground state manifold is closed 
and completely connected - no disconnected subsets of 
the manifold exist. 

We now assume that these 8 fields are slowly varying 
with position on the lattice, so we can treat them us- 
ing a continuum limit. However, these 8 fields will still 
be required to respect the symmetry of the underlying 
lattice. 

The 8 slowly varying fields 4> u . a transform under a par- 
ticular 8-dimensional irreducible representation (irrep) of 
the PSG. In fact, we prove in appendix IF1 that this is the 
minimum dimension of a representation of the PSG. It is 
convenient to change the field basis to 





C0 + C1 




C0-C1 


Ou = 


V2 ' 


<f>2u = 






C2 + C3 




6 -6 


lu = 


V2 ' 


<t>3u = 






6 + 6 




6-6 


>0d = 




<j>2d = 


V2 




6 + 6 




6 - 6 


>ld = 


V2 ' 


<j>3d = 


v/2 



(38) 



This basis realizes a "permutative representation" of the 
PSG, in the nomenclature of Refi That is, the sym- 
metry operations of the PSG act on these fields by a 
combination of permutations and simple diagonal phase 
rotations. See Eqs. I|D1|) - (|D6(I . Because of this structure, 
the action takes a particularly simple form in this basis. 



C. Effective low energy action 

Our goal is to describe a low-temperature condensate 
phase of the monopole defects. Limiting the discussion to 
zero temperature, we set out to construct a phenomeno- 
logical Landau- Ginzburg (LG) action to access the con- 
densate phase. We concentrate on the ground state man- 
ifold of the monopole defects, ignoring any higher energy 
modes, and construct a low energy effective continuum 
action in the 8 field components of l|37|) . The various 
terms allowed in the action must be invariant under the 
PSG in the 8-dimensional representation (|D1|) - (|D6|I . As 



11 



the effective action will live in a 3+1 dimensional space- 
time, we seek terms only up to quartic order in the field 
operators - any higher order terms will be irrelevant in 
the renormalization group sense. 

To quadratic order only one invariant exists, 



,. 2 ) , 



(39) 



i\ 2 > 



(40) 



a typical mass term, where the indices enumerate j = 
0,1,2,3. 

At quartic order we find 4 invariants 

(Ox) 2 = (|cl 2 + |el 2 ) 2 , 
e 2 = £lOI 2 IOI 2 + ElOI 2 l£ 

e 3 = |cTlH 2 ) 
e 4 E/r,/ w C;.;"/Cs , 

ijkl 

where the sums over indices are always over 0,1,2,3 un- 
less otherwise stated, and the vector notation implies a 
4-component vector. We shall not specify the k tensor 
explicitly, as it is long and complicated. We leave it to 
Appendix [E] 

The new set of independent terms obeys numerous con- 
tinuous symmetries in a rather transparent manner (for 
brevity we do not specify the discrete symmetries of these 
terms): 

• Oi is invariant under a full U(8) symmetry, as it is 
just the norm of an 8 dimensional complex coordi- 
nate vector. 

• G>2 is invariant under a [t/(l)] s symmetry, since it 
depends only upon the magnitudes of the field com- 
ponents in the new basis. We may change freely 
and independently the phase of each field compo- 
nent. 

• O3 is invariant under a £/(4) x C/(4) symmetry, be- 
cause it involves only the norm of two 4 dimensional 
complex vectors. 

• O4 obeys a [£/(l)] 2 symmetry group corresponding 
to the following transformation 



Vj 








(41) 



The action as a whole is invariant under a [£/(l)] 2 sym- 
metry, governed by the 64 invariant. 

The only microscopic symmetry of the action is the 
dual U(l) gauge invariance, i.e. identical phase rotations 
of all field components. This changes just the overall 
phase of the monopole wavefunction, and has no physical 
consequence (it is a true gauge). The other "staggered" 
U(l) rotation (with 6 = —A in Eq. 1411) ) is accidental, 



occurring only in the proximity of the critical point at 
which truncation to the quartic action is a good first ap- 
proximation. It will be broken if we include sufficiently 
high order terms in the action. Our investigation con- 
cluded that the staggered U(l) symmetry persists at 6th 
order, but ultimately breaks down to a discrete symme- 
try (Z4) at 8th order. The remaining discrete symmetry 
Z4 can be identified as part of the lattice PSG. 

Finally, the most general low energy effective action up 
to quartic order is 



S 



•7i 



{\(d^-ia^)Cj\ 2 + \(d fi -ia fl )^\ 2 + ^ 

e 2 + 72 e 2 + 73 e 3 + 7464 + ^f 2 } , 



(42) 



where a sum over j, as well as integration over spacetime 
are implied, and -^iF 2 is the Maxwell term in the action, 
with F^y = d^a^—d^a^. Here a is the continuum version 
of the vector potential a appearing in Eq. J57J. The 
7j are phenomenological couplings undetermined in our 
theory. 

Higher order terms can be ignored for a renormaliza- 
tion group (RG) analysis (in 3 + 1 dimensions, they are 
irrelevant in the RG sense, and can be treated perturba- 
tively when necessary), but must be taken into account 
in a mean field analysis to which we now turn. 



V. MEAN FIELD THEORY 

We now turn to analyze our effective action (|42|l . Using 
mean-field theory (MFT) we can find the various phases 
of this action. Replacing the fields with their average 
values, we obtain a Lagrangian density C that must be 
minimized with respect to the field averages. 

At first we ignore the 04 term, since it has the lowest 
symmetry. We shall examine only the ordered phases 
which occur for weak 74 (and higher order coefficients 
75 ! 76 when necessary - see below) for simplicity. Scaling 
all the couplings to that of Q\ 2 we get 

M 2 = - — , 
7i 



Pl = 


72 




7i ' 


02 = 


^ + 2, 




7i 


03 = 


74 




7i ' 



(43) 



The mean-field Lagrangian density with (3% = 74 = is 

£/7i = ^-m 2 ic? + (id 2 ) 2 +/?iE ic/Kii 2 j + (c - 

+ /3 2 |Cl 2 |£?- 

(44) 
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FIG. 8: (Color online) Mean field theory phase diagram of 
simplified free energy, fiz = 0. The phases exist only within 
the bounds of the dashed lines. Different phases are delimited 
by continuous lines. 



We first find the mean-field minima with M 2 > in the 
01,2 plane. The resulting phase diagram is illustrated in 
Fig. |S| We label the "phases" by the number of non-zero 
components of £ and C, ordered from largest to smallest, 
with the notation N\ + iVjj. The "phases" are: 

• X > , 2 > 2, 

In this region, only one field component condenses, 
so this region is labeled 1 + 0. 



01 > , -2 < 2 < 2, 

Here one field component of each set 

denses, so this is a 1 + 1 region. 



con- 



• Pi<O,02 > 2(1 + ^1), 
Here all the field components of one of the 2 sets 
condense, denoted 4 + 0. 

. X > 0,-2(1 + /3xl) <0 2 < 2(l + ft|), 
Here all the 8 field components are non-zero, so it 
is a 4 + 4 region. 

It should be pointed out that specifying the number 
of non-zero components of the complex fields does not 
necessarily specify a unique charge-ordered phase of the 
problem. Physically, these charge-ordered states break 
only discrete space group symmetries, and so are ex- 
pected to lead to only a discrete degeneracy of ground 
states. In our dual gauge formulation, each physically 
distinct state is represented by a "cycle" in the monopole 
field space, along which the phase of all the fields is varied 
together by a uniform angle which can go from to 27T. 
This variation has no physical significance since this is 
a gauge symmetry. Other continuous phase freedoms of 
the mean field solutions, which do not vary the phases of 
all components identically, are however not a pure gauge 
freedom. Such a (artificial) freedom can only be a result 
of truncating the dual action to quartic order. 

In fact, only the 1 + phase lacks any such "emer- 
gent" phase freedom. In particular, the different 1 + 



states obtained by choosing each of the 8 possible com- 
ponents non-zero are physically distinct, and the phase 
of this non- vanishing field is pure gauge, with no physical 
significance. Thus, the 1 + states are eight-fold degen- 
erate. The 03 term vanishes for such states and does not 
change this conclusion. These states exhibit an enlarged 
unit cell, comprising 4 unit cells of the underlying lattice, 
and containing 16 pyrochlore lattice sites (In fact, we will 
find that all our mean-field states have an enlarged unit 
cell. The amount by which the unit cell is enlarged can be 
computed by explicitly checking how a given mean-field 
state transforms under an arbitary translation using the 
PSG.) The Bravais lattice formed by these enlarged unit 
cells is simple cubic (SC). It is noteworthy that this phase 
is the R state found in Ref. Ill 

The other 1 + 1, 4+0, and 4+4 solutions retain at least 
some physically meaningful phase freedom even when 03 
is included in the action. Consider first the 4+0 case. For 
such configurations, the 03 term vanishes, so the phases 
of all four non-zero fields remain free at quartic order. 
Only the overall phase of the four fields is gauge, and 
hence will remain undetermined - and unphysical - to 
all orders. To fix the remaining three phases, we must 
consider higher order terms. In this case, 6 th order is 
sufficient, and there is one term which resolves the con- 
tinuous degeneracy: 

C 6 = 75Re{ e -" /4 [(Co) 3 CiC 2 C3-(Ci*) 3 CoC 2 C3 

-(c 2 *) 3 c CiC 3 + (a) 3 c CiC 2 ] - (c -> D}- (45) 

In the N = 4+0 phase for 75 > 0, one minimum energy 
state is 



l,e +l T- ; e H 



(46) 



It is part of a set of 792 total states, which can be gener- 
ated from this one by the PSG. For 75 < 0, one configu- 
ration is 



c = (i 



M) 



(47) 



which is part of a set of 384 degenerate states connected 
by the PSG. The enlarged unit cell comprises 16 unit cells 
of the original lattice, and contains 64 pyrochlore lattice 
sites. The "supercells" are arranged in a body-centered 
cubic (BCC) Bravais lattice. 

In the 1 + 1 case, there is a single undetermined phys- 
ical phase, an opposite phase rotation of £ and £. This 
is a direct result of the "staggered" U(l) symmetry of 
the action at 6 th order. Any such solution breaks this 
staggered U(l), and so making such a rotation gives a 
different solution. The staggered U(l) symmetry is bro- 
ken only at 8 th order, by 8 distinct terms. In the 1 + 1 
case, only one of these terms is non- vanishing. It takes 
the form 



£ 8 = 



76 Re{[(C *) 4 -(Cr) 4 -(C2) 4 + (C3) 4 ] 



4 4 ]}. 



(48) 
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The nature of the phase depends upon the sign of 76. In 
each case, there are 64 physically distinct ground states: 
4 choices each for the non-zero component of and 4 
inequivalent phase minima of Eq. (|48J) . 

Of the many states we mention one configuration for 
76 < 



Co = £0 = 1 , 
and one configuration for 75 > 



Co 



r 5o 



1 



(49) 



(50) 



The enlarged unit cell comprises 16 unit cells of the orig- 
inal lattice, and contains 64 pyrcohlore lattice sites. The 
"supercells" are arranged in a simple hexagonal Bravais 
lattice, with a ratio of - = 4VL5 in standard notation. 

In the 4 + 4 case the relative phases in each quartet are 
determined at 4 th order by /?3, and the same "staggered" 
£7(1) symmetry as in the 1 + 1 case remains. 

For 74 > one possible configuration is 



C = (l,e+^,e +l T5,e-^) , 

£ = e +lXi (l,e +4 T5,e+^,e-^) , 

and for 74 < a possible configuration is 



(51) 



C= U,e 



1, e 1 4 , e 1 4 , e" 



(52) 



Many (over 10000 in each case) other configurations are 
possible, and connected via the PSG to this state. 

The "staggered" U(l) symmetry is lifted by various 8 th 
order terms. Depending on the various couplings of the 
8 th order terms, we get that the "staggered" U(l) breaks 
into the same two Z4 multiplets as the 1 + 1 state does. 
We shall not belabor the details of all 8 of these terms. 
(|48|l alone suffices to access both multiplets. For 76 > 



A4 — An 



7T 

4 



and for 76 < 



A4 — An = n— . 



(53) 



(54) 



Finally, in any one of these cases, the enlarged unit cell 
comprises 4 x 4 x 4 — 64 unit cells of the original lattice, 
and therefore contains 256 pyrochlore lattice sites. The 
"supercells" are arranged in a face-centered cubic (FCC) 
Bravais lattice, with the primitive lattice vectors taken 4 
times larger relative to the original FCC. 

Using a direct mapping between the monopole defect 
density and the spin density variations described in the 
next section, we can depict the latter in Figs. II 11191 



VI. ORDERING PATTERNS 

Having found the various allowed monopole defect con- 
densate phases in the abstract order parameter space by 
algebraic considerations, we want to identify the physical 
spin (Sf) ordering patterns in each case. 

The density of the monopole defects is a varying scalar 
density. As it is the only spatially varying scalar in the 
problem, the spin density must have similar spatial vari- 
ation. More precisely, both densities must obey the same 
symmetries, 



(55) 



where p(f) is the monopole defect density. The monopole 
defect density can be found by taking the wavefunction 
(|3T|l squared, 



p(f) = |*(r)| 5 



(56) 



The resulting density is now a function of the field val- 
ues 4> u ,a- By inserting the values for the different MFT 
phases, we can recover the density. 

In some phases, however, the symmetry breaking is not 
manifest in the scalar density, but rather in the current 
or kinetic energy: 



J r r , = ih*{f)^{f) e - ia ^' - c.c) 
K r ,T> = (^*{r')^{r)e- ia -y + c.c) 



(57) 



Both the current density and the local kinetic energy 
can be encoded in a complex valued vector, 



v rr i = ^*(r)^(r)e~ 



(58) 



The imaginary and real parts will give us (half) the cur- 
rent density and the local kinetic energy, respectively. 

Each plaquette in the dual diamond lattice corresponds 
to a pyrochlore lattice site at the center of the plaquette. 
Any monopole defect "object" we can define on the dual 
plaquettes is also defined on the direct pyrochlore lat- 
tice sites, and encodes the symmetry of the MFT phase. 
Therefore, the function must be "similar" to the spin 
density on these sites, in the sense of giving the correct 
symmetry of the latter. An appropriate function is a 
loop integral (curl) of the complex current around the 
plaquette 



vppi ~ no = Tii 



(59) 



ff'GO 



We can formalize this argument by considering 
Maxwell's equations, with magnetic monopoles 



-* dB -> 
curli? = — — + J b , 
at 

— dE -> 
curlB = +— + J e , 
at 

div E = p e , 
div-B = pt, . 



(60) 
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FIG. 9: (Color online.) Notation for the spin density pictures. 
Triangles with a site at the center represent up-pointing tetra- 
hedra. Triangles with no site at the center represent down- 
pointing tetrahedra. 



where the magnetic monopole density and current are 
denoted with a subscript b. In a static system integrating 
the first equation over some surface we get, by Stokes 
theorem, 



J b -d£ = 



cm\E-d£= / E-dA. (61) 
c Js 



From this last expression it is evident that the loop in- 
tegral of the monopole defect current gives the electric 
flux through that loop. The lattice version of the electric 
flux is a constant plaquette "area" times the electric field 
penetrating perpendicular to that plaquette 



Jff' — E a b 

rr' E O 



(62) 



In conclusion the variations in the spin density can be 
related to the loop sums of the monopole defect current 
around plaquettes of the dual diamond lattice. Armed 
with this knowledge we can plot a function to show the 
spin density variations. We plot the "spin density" so de- 
fined for each of the phases obtained in mean field theory 

in Figs, mm 



VII. RG ANALYSIS 

In this section, we briefly consider the effect of fluctu- 
ations on the mean-field critical behavior of our effective 
action. Our primary focus in this paper is not on quan- 
tum critical phenomena, but rather on the nature of the 
ordered phases close to the U(l) spin liquid state. These 
results for the ordered phases are independent of the con- 
tents of this section. 

By simple power-counting, the problem of a general- 
ized "Ginzburg-Landau" theory in 3+1 dimensions (with 
a many-component "superconducting" field ife) is in its 
upper critical dimension, so one expects either of two pos- 
sibilities. One possibility is that the Gaussian fixed point 
is marginally stable, and mean-field behavior is correct 
up to logarithmic factors. The other possibility is that 
the Gaussian fixed point is marginally unstable, and the 
true critical behavior is a strong coupling problem; most 
probably, such flows to strong coupling indicate a weak 
fluctuation-induced first order transition. 

Here we follow Balents et aii£, which generalized the 
calculations of Halperin, Lubensky and MsM> and Brezin 




FIG. 10: (Color online.) The spin density variations in the 
1 + phase - a 3D image. 1 field component only has a non- 
vanishing expectation value. Here only minority sites in one 
kagome plane and in the triangular plane "above" the page 
are shown. Those in the triangular plane "below" the page 
are omitted to keep the image uncluttered. 




FIG. 11: (Color online.) The spin density variations in the 
1 + phase - same phase as in Fig. HUI for comparison. This 
phase has an enlarged unit cell of 2x2x1 = 4, in a simple 
cubic Bravais lattice. 



et aii«a, and consider a general g-component [[/(!)] action 
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FIG. 12: (Color online.) The spin density variations in the 
1 + 1 phase, for 75 < 0. One each of the four £ and £ fields has 
a non-vanishing expectation value. The expectation values 
have identical magnitude. This phase has an enlarged unit 
cell of 2 x 4 x 2 = 16, in an simple hexagonal Bravais lattice. 



FIG. 14: (Color online.) The spin density variations in the 
4 + phase, for 74 > 0. The 4 field components of either the 
£ or £ set have a non-vanishing expectation value of identical 
magnitude. This phase has an enlarged unit cell of 16, in a 
BCC Bravais lattice. 



FIG. 13: (Color online.) The spin density variations in the 
1 + 1 phase, for 75 > 0. One each of the four £ and £ fields has 
a non-vanishing expectation value. The expectation values 
have identical magnitude. This phase has an enlarged unit 
cell of 2 x 4 x 2 = 16, in an simple hexagonal Bravais lattice. 



9-1 



So = [ d D r{ ]T \(d, - ia»)w\ 2 + s\w\ 2 + 7^F 2 } , 

J I) n 



e=o 



i r q ~ x 

£,m,n,i=0 



(63) 



Here we have written the theory for a general space-time 
dimensionality D. For the quantum critical point of in- 
terest, D = 3 + 1 = 4 total space-time dimensions. For 
this very general action the RG flows obtained by an e- 



FIG. 15: (Color online.) The spin density variations in the 
4 + phase, for 74 > 0. The 4 field components of either the 
£ or £ set have a non-vanishing expectation value of identical 
magnitude. This phase has an enlarged unit cell of 16, in a 
BCC Bravais lattice. 



expansion arc 
de 2 2 Cq 



di 
dU£ r , 



d( 



= (e + 2rj)ui m . nk - a 2 e 4 C [6i n S mk + 6e k Sr, 



ij 

2 



r] = die C 



(64) 



with spatial dimension d — 3 — e. Here r\ = aie 2 C is the 
anomalous dimension of the fields. The constants C, 0,1,2 
were calculated^ in an e-expansion. For our case e = 
since our model is in 3 + 1 spacetime dimensions. Thus 
the e-expansion results hold exactly: C — g^-, a\ = 
3, a 2 = 6. 
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FIG. 16: (Color online.) The spin density variations in the 
4 + 4 phase, for 74 > 0, 76 > 0. All 8 field components have a 
non-vanishing expectation value of identical magnitude. This 
phase has an enlarged unit cell of 4 x 4 x 4 = 64, in a FCC 
Bravais lattice. 



FIG. 18: (Color online.) The spin density variations in the 
4 + 4 phase, for 74 > 0, 76 < 0. All 8 field components have a 
non-vanishing expectation value of identical magnitude. This 
phase has an enlarged unit cell of 4 x 4 x 4 = 64, in a FCC 
Bravais lattice. 



FIG. 17: (Color online.) The spin density variations in the 
4 + 4 phase, for 74 < 0, 76 > 0. All 8 field components have a 
non-vanishing expectation value of identical magnitude. This 
phase has an enlarged unit cell of 4 x 4 x 4 = 64, in a FCC 
Bravais lattice. 



FIG. 19: (Color online.) The spin density variations in the 
4 + 4 phase, for 74 < 0, 76 < 0. All 8 field components have a 
non-vanishing expectation value of identical magnitude. This 
phase has an enlarged unit cell of 4 x 4 x 4 = 64, in a FCC 
Bravais lattice. 



We take directly ip — > £, £ in the action 1)42(1 to make 
contact with the couplings 7^ 

^mm;mm — 47l ; 

U m n;mn = 4(72 + ~7l) II / III , 

2 (65) 

Umn-mn = (73 , 

where the indices € 0...3 and the over lined indices de- 
note the £j quartet of field components. All other Uijki 
couplings vanish in our special case. This "structure" 
encodes in it the specific symmetry group of our lattice 
model. Therefore, the RG flow equations should preserve 
this general structure - only 4 independent couplings, 
and the K-tensor must retain its structure. This provides 
a good check that we have indeed included all possible 



quartic terms allowed by the symmetry of the problem. 

Some manipulation gives us the flow equations for our 
4 couplings: 

^ = _ C ( 3e 4_ 6e 2 7l+87l 2 + 

at x 

16 71 72 + 16 72 2 + 8 7i 73 + 2 7 3 2 ) , 

^=C (6e 2 7 2 - 2 7 i 7 2 + 272 2 - 74 2 ) , 

f (66) 

^=2C(l6 7 2 2 + 3 e 2 73- 

6 71 73 - 8 72 73 + 73 2 - 3 7 4 2 ) , 

^ =2C (3e 2 - 671-273-74) 74 ■ 

The only fixed point allowed by the RG equations is a 
trivial fixed point with all coupling strengths vanishing 



7i = 0, e = 0. 

The stability or lack thereof of such coupled non-linear 
differential equations is not obvious. We would like to 
know if there is a subspace of codimension zero of the 
four dimensional phase space (we may project it onto 
the e = plane since the evolution e > is clearly 
monotonically decreasing toward zero) in which all cou- 
plings scale toward zero. While we have not been able 
to prove this is not the case, our numerical and ana- 
lytical investigations have found no such stable regime. 
We did find a few specific fine-tuned stable solutions, 
but these were all unstable to infinitesimal perturbations. 
Thus we believe that the mean-field critical behavior is 
always destabilized by fluctuations. We expect this prob- 
ably signals a fluctuation-driven weakly first order tran- 
sition. Intuitively, this is a result of the gauge field fluc- 
tuations^, which lead to attractive interactions amongst 
the monopoles, driving bound state formation. Note that 
this result is reliable only for the zero temperature tran- 
sition with D = 3 + 1 = 4, the upper critical dimension, 
where the perturbative renormalization group treatment 
is justified. In D = 3 + 0, the same field theory may 
well have non-trivial stable (if one parameter is tuned to 
criticality) fixed points. We will return to this point in 
the discussion below. 
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FIG. 20: (Color online.) Schematic phase diagram of the 
spin-dimer model in the v-T plane. At zero temperature, we 
suppose there is a direct transition from a U(l) spin liquid to 
the R state at v = v c < 1. The thick (red online) portion of 
the T > phase boundary is first order, while the thin (blue 
online) boundary denotes a possible non-trivial "non-LGW" 
second order transition at higher temperature. 



VIII. DISCUSSION 

In the preceding sections, we have presented a system- 
atic study of the zero temperature ordered phases proxi- 
mate to a "half-polarized" (from the 3:1 constraint) U(l) 
spin liquid on the pyrochlore lattice, based on a projective 
symmetry group (PSG) analysis of the monopole excita- 
tions of the liquid state. One of these states would be 
expected to occur on reducing v in the spin/dimer model 
starting from a value (< 1) within the spin liquid state. 
The ordered phases, determined at the mean-field level, 
all break discrete symmetries of the pyrochlore lattice, 
and in particular we find a generic unit cell enlargement, 
with the minimal cell size 4 times that of the underly- 
ing pyrochlore lattice, and a maximal unit cell 64 times 
larger. The simplest "R" state, with the smallest unit 
cell, is the same one which was found to be the ground 
state in the "classical" spin/dimer limit v — ► — oo. It 
therefore seems likely that the spin/dimer model exhibits 
only two phases for v < 1, and hence a direct quantum 
phase transition from the U(l) spin liquid to the R state. 
If this is indeed the case, one may contemplate the pos- 
sibility that HgCr204 (or other pyrochlore chromates) 
might be close to this quantum critical point. 

While previous worki^ and the main text has focused 
almost entirely upon the ground state properties of the 
spin/dimer model, it is interesting to consider the more 
general problem at non-zero temperature T > 0. A 
schematic phase diagram in the v-T plane is illustrated 
in Fig. |201 We will focus on the region v < 1, and will not 
discuss the physics of the thermal ensemble of "frozen" 



states occuring for v > 1. First, we note that the R state 
has a gap to all excitations, and breaks only discrete lat- 
tice symmetries. We therefore expect that the R state 
will persist at T > up to some non-zero critical tem- 
perature T c (v) > for v < v c . By contrast, the U(l) 
quantum spin liquid ground state breaks no symmetries. 
Therefore, in the dimer model, no transition is expected 
for v c < v < 1 as T is increased from zero to infinity. We 
thus expect just the single phase boundary emanating 
from the quantum critical point, shown in Fig. I2UI 

It might appear from these observations that the ther- 
mal ordering transition at T c (v) should be a rather or- 
dinary one, described by the usual Landau-Ginzburg- 
Wilson (LGW) approach based on the order parameter 
of the R state. In fact this is incorrect, and T > prob- 
lem is rather more interesting. To see this, let us con- 
sider the "paramagnetic" state obtained for very high 
temperature in the dimer model. The physics of such 
classical, infinite temperature dimer models (and other 
similarly constrained models) have been considered by 
several authora&SSiSi. As shown in these works, due to 
the dimer constraint, even at infinite temperature the 
dimer model has non-trivial power-law "dipolar" corre- 
lations. Such dipolar correlations are not captured by the 
conventional LGW theory which retains only the order 
parameter. 

To understand these dipolar correlations, and the 
a proper formulation of the phase transition in the 
spin/dimer model, it is instructive to return to the 
monopole field theory in Ea. (|63|l . Let us rewrite this 
effective action making space and time coordinates sep- 
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arate and explicit: 
S = 

r r c q ~ 1 

/ d * r d i 5Z ~ *«oM 2 + |(V - ia)^| 2 + S |^| 2 

+ i(^-S 2+ 2?P x5 )T 

1 f f 13 9-1 

•'° e,m,n,i=0 

Here /3 = (fc B T) _1 is the inverse temperature. To derive 
a theory of the thermal phase transition at T > 0, we use 
standard logic to proceed from Eq. (|67|l . In particular, 
at T > 0, imposing periodic boundary conditions (as 
demanded by the trace defining the quantum statistical 
mechanical partition function) in imaginary time leads 
to a set of discrete bosonic Matsubara frequencies u n = 
27rn//3, with integer n. Because of the time-derivative 
term, all modes with uu n ^ have enhanced "masses" 
relative to the zero Matsubara frequency mode, and can 
be integrated out. Practically speaking, this amounts to 
assuming the order parameter is constant in imaginary 
time, d T (f£ = 0. Similarly, we may take d T a = 0, and by 
a choice of gauge, ao — 0. Carrying out this procedure, 
we find Sq + Si = T '/k B T, where T = + T\ is an 
effective classical free energy: 

r < q ~ x 

r Q = / d 3 r[J2\(V-ia)<p £ \ 2 + s\<p e \ 2 
^ 1=0 

+ i(^ X<3 ) 2 }< 
If 9-1 

•T 7 ! = 4 / H Uem;ni l Pt l Pm l Pn ( Pi ■ ( 68 ) 

Eq. I)68|l is precisely the classical "Ginzburg-Landau" 
free energy for a multi-component superconductor in 
three dimensions, with the quartic interaction ue m ;ni de- 
termined by the PSG. The interpretation of this result is 
quite simple. In the U(l) spin liquid, the monopole is a 
well-defined, bosonic particle excitation, and carries the 
"magnetic" gauge charge. The magnetic gauge charge is 
conserved in the theory. As already discussed, the or- 
dered (R) state at zero temperature is understood as a 
condensate of these monopoles. In fact, within the quan- 
tum dimer model, the bosonic monopole can condense 
at a non-zero temperature, just as in an ordinary Bose- 
Einstein condensate. Thus Eq. (|(j8fl is nothing but the 
"classical" free energy describing the "superfluid" tran- 
sition of this monopole. Because it carries a non-zero 
magnetic gauge charge, it is coupled to the electric vec- 
tor potential a, just as the electric charged Cooper pair 
condensate is coupled to the magnetic vector potential in 
ordinary superconductivity. 

Thus we are led to this remarkable and unconventional 
description of the T > phase transition between the 



paramagnet and the R state. Without the 3:1 spin/dimer 
constraint, this transition would certainly be expected to 
be governed by conventional Landau theory. What is the 
nature of the unconventional transition? According to 
the d = 4 — e expansion approach (discussed in Sec. IVIIfl . 
this transition is fluctuation-driven first order. However, 
this conclusion is known to be often incorrect for the 
physical case of Ginzburg-Landau transitions in three di- 
mensions. In particular, a class of related models, N- 
component Ginzburg-Landau theories with U(N) sym- 
metry, have been investigated in a number of cases. First, 
for sufficiently large N, these transitions can be shown 
to be second order in an expansion around N = oo. Sec- 
ond, for N = 1, a duality transformation has been used 
to demonstrate that the transition can be continuous, in 
the inverted XY universality class^. A similar duality 
analysis, in conjunction with numerics, has been used 
in arguing for continuous critical behavior for N — 2 34 . 
It therefore appears quite likely that continuous critical 
behavior is possible in these models for any N. By anal- 
ogy, continuous critical behavior of the theory of Eq. 
seems quite possible. We thus suggest that the param- 
agnetic to R state transition in the spin/dimer model 
constitutes a novel non-LGW universality class. Note 
that, since the RG analysis in Sec. I V lll concluded that the 
T = quantum phase transition is weakly first order, the 
phase boundary for very small but non-zero temperature 
must remain discontinuous. We therefore expect a mul- 
ticritical point separating on the T > phase boundary 
separating first order from continuous non-LGW critical 
behavior. The continuous and first order portions of the 
phase boundary are shown in Fig. l20l bv thin (blue online) 
and thick (red online) lines, respectively. It would be of 
considerable interest to investigate this classical phase 
transition in the the spin/dimer model numerically. This 
could be done on a purely classical dimer model, with 
K = 0, and so would require only classical Monte Carlo 
methods. 

As remarked above, without the 3:1 spin/dimer con- 
straint, the classical phase transition would certainly be 
expected to be described by LGW theory. For any mi- 
croscopic model, such as the spin-3/2 Heisenberg model 
of Reflla. the 3:1 constraint is not expected to be ex- 
actly obeyed. At zero temperature, however, pyrochlore 
tetrahedra violating the 3:1 constraint are gapped exci- 
tations, and are not present in the ground state on the 
plateau. Consider configurations in which a single tetra- 
hedron violates the constraint by having either zero or 
two minority spins instead of one. In the mapping to the 
gauge theory, these particular excitations can be viewed 
as states with "electric" gauge charge ±1 (relative to the 
static background gauge charge of the plateau states) on 
the tetrahedron in question. They also carry physical 
spin AS Z = ±3/2 relative to the plateau states. This 
follows because the total spin can be written as 

t 
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the factor of 1/2 being required since each spin is con- 
tained in two tetrahedra. These excited states can 
therefore be viewed as fractional spin excitations, or 
"spinons" , in this particular example. In any case, be- 
cause they cost only finite energy, there will be a non- 
vanishing concentration of such electric gauge charges at 
T > 0, due to thermal activation. The typical separa- 
tion of the electric gauge charges is expected to behave 
exponentially at low temperature, A ~ exp(A/fc B T), if 
A is the gap to the electric charged particles. This 
has several effects. First, the dipolar correlations of 
the paramagnetic phase of the spin/dimer model will 
cross over to the usual exponential ones of an ordinary 
paramagnet, for lengths larger than A. Moreover, the 
"plasma" of electric charges is expected to give rise to a 
linear confining potential between oppositely (magneti- 
cally) charged monopoles. This will bind the monopoles 
and anti-monopoles into gauge-neutral pairs, the radius 
of this bound state being at least as large as A. It is 
only these pairs that can Bose condense. When A is 
large, there is a cross-over behavior. On approaching 
the paramagnetic to R state transition, the correlation 
length grows in a manner first governed by the non-LGW 
theory of Eq. I|68|) . Once the correlation length exceeds 
the monopole-antimonopole binding length, these pairs 
may be considered well-formed, and the critical behav- 
ior changes to that described by Bose condensation of 
the pairs. The creation/annihilation operators for the 
monopole-antimonopole bound states is expected to be 
proportional to 



<Pt<Pt' 



(70) 



Since Lp\ and ipi transform under conjugate representa- 
tions of the PSG, and hence is gauge-neutral, it 
transforms not under the PSG but simply the ordinary 
lattice space group. Hence ^u 1 can be decomposed into 
irreducible representations of the space group, which are 
precisely the usual Landau order parameters! Thus the 
critical behavior sufficiently close to the phase bound- 
ary, when the 3:1 spin/dimer constraint is not rigidly 
enforced, is indeed governed by LGW theory as expected 
on general grounds. The conventional LGW analysis is 
straightforward, and will be presented in Ref.f3^ It pre- 
dicts that the paramagnetic to R state transition should 
be first order, due to the presence of a cubic invariant. 
For k B T <C A, it will be weakly so along the thin (blue) 
portion of the phase boundary in Fig.[3U| Of course there 
is no such crossover in the spin/dimer model, for which 
the constraint is rigidly enforced. 

In summary, we have studied the phase structure of 
a spin/dimer model on the pyrochlore/diamond lattice. 
It contains an interesting quantum paramagnetic "U(l) 
spin liquid" phase at zero temperature, as well as an or- 
dered "R" state at T > 0. We derived a novel monopole 
field theory to describe the quantum and classical phase 
transitions between the quantum and classical paramag- 
nets and the R state. Prior work in ReflT3l indicates that 
the half-polarized magnetization plateau of the spin-3/2 



Heisenberg antiferromagnet on the pyrochlore lattice is 
described by this model, with a coupling constant that 
places it in the R state at low temperature. It is in- 
teresting to contemplate the possibility that the experi- 
mental materials CdCr204 , HgCi'204 might be near the 
quantum phase transition of the paper to the spin liquid 
state. It would indeed be exciting were some homolo- 
gous material to actually realize the U(l) spin liquid in 
its plateau ground state. To this end, we note that, were 
a spin- 1/2 pyrochlore antiferromagnet to be realized ex- 
perimentally, the quantum effects would be significantly 
further increased, and a U(l) spin liquid might well be 
expected theoretically. We leave such tantalizing possi- 
bilities as open questions. 
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APPENDIX A: GROUND STATE EIGENMODES 

The diamond lattice is bipartite and can be thought of 
as being made of 2 FCC sub-lattices one shifted from the 
other by the vector b = — \(a>\ +02 + 03) where Sj are 
the FCC primitive vectors. We introduce the notation 



771 (n) = tp(f) , 

772 (n) = ip(r + b) 



(Al) 



where r denotes sites on an FCC lattice, r + 6 are the 
sites of the second FCC lattice comprising the diamond 
lattice, and ft are the "index" space coordinates on the 
FCC lattice. 

Defining the Fourier transform of the wavefunction 77, 
and its inverse as 



r\ (ft) = 



d 3 k 
(2+p 



f)(k) ■ e 



r](k) = 



keBZ 



(A2) 



V(n) , 



we can write out the eigenstates in a compact manner 



iu=0...3 
a=l,2 



(A3) 



fi=0...3 
a=l,2 



where |0) denotes the vacuum state, and 
Pn,v =p + f-iQ + -ue , 
= d-p + fiQ + -ve , 



(A4) 
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where p — (5, j, f ) and e? = (5, 0,7r) . The coefficients 
are 



(A5) 



where c^ ' = cJ^q are the coefficients for the v — Q state. 
Finally, the un-normalized coefficients cjf^ are: 

■i 1 



1,-*, 



c< 2 > = 



1 + v/2' l + \/2/ ' 
3 - i l-2i 1 + 3* 2 + i 



lo'v^ + yio' \/T0(l + V2)' V5 



(A6) 



APPENDIX B: PROJECTIVE SYMMETRY 
GROUP (PSG) 

We start by enumerating some diamond lattice symme- 
tries. Since the diamond lattice is made up of two FCC 
lattices, it inherits the translations of the FCC lattice 

ij : r — ► f + 3j . (Bl) 

Next we consider a 3-fold rotation symmetry 

?3 : {x,y,z) — > (z,x,y) , (B2) 

and a 2-fold rotation 

r 2 : (x, y, z) — ► (x, -y, -z) . (B3) 

There is also a reflection symmetry 

i : (x,y,z) — ► (x,z,y) . (B4) 

Finally we have an inversion symmetry that effectively 
swaps between the two FCC sublattices 



p : r 



(B5) 



To construct the PSG we attach gauge transformations 
to each symmetry operation. This is most easily done in 
the "index" space coordinates, and we now list the action 
of the PSG operations on the wavefunctions. First are 
the three lattice translations 



fj : n(f£) — > f](n + a 3 )e- lA ^ H) 

where 

A j(n) =-(£*• n)Qj . 
The inversion symmetry 

V : T)(n) — ► a x ■ r)(-n)e- lKc ^ 
where a x is the x Pauli matrix, and 



A c (n) 



-d-n 



(l,0,-2)-n 



(B6) 
(B7) 
(B8) 
(B9) 



The 3-fold rotation 

K {3) : 77(7?) — ► V (H [3) -n)e~^ n \ 
with the gauge transformation 

A7J3 (ft) = ft ■ Bi ■ ft + S ■ ft , 

where 

A 1 2\ 

6 = -t (1,2,1) 



iA^(n) 




Bi-- 

The 2-fold rotation 



(BIO) 



(Bll) 



(B12) 



K 



(2) 



772(7?) > m(R( 2 ) -n + di)e 



-iA^(n) 



(B13) 



The gauge transformations for the 2 spinor components 
are different 



where 



A^ (77) = 77 • B 2 ■ n + Xi • 77 , 
A y K ' 2 (ft) = ft ■ B 2 ■ ft + A 2 • 77 , 



B 2 = 




(B14a) 
(B14b) 



(B15) 



and 



Ai 



\ (-1,1, 2) A 2 = A 1 + |(1,2,1) . (B16) 



Finally, the reflection symmetry 

X: 77(77) — ► tf(I -n)e- ih * (K) 
and the gauge transformation is 

Aj(77) = 77 • B 3 ■ ft + £ • 77 , 

where 



B 3 = - 




(1,1,1) 



(B17) 



(B18) 



(B19) 



One can take the Fourier transform of these symme- 
try operations, to find their action on r)(k) and then de- 
duce how the 8 ground state eigenmodes transform under 
them. The transformation rules I|D1I) - (|D6(I are the result 
of that analysis. It is a good check that the ground state 
manifold is invariant under every one of these symme- 
tries. Although the ground state manifold need not be 
completely connected by the symmetry operations, in our 
case it is. 
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APPENDIX C: GROUP ALGEBRA 

The diamond lattice symmetry group, described by the 
reduced set of operations introduced in the text, obeys 
the algebra 





r 3 3 


— 


1 






r 2 2 


= 


1 






P 2 


= 


1 






i 2 


= 


1 




r 3 ■ P 


■ r 3 ~ 1 ■ p 


= 


1 




% ■ L ■ 

1 3 


f.- 1 . f.-l 

i "J 




1 






(p ■ hf 




1 




J 3 


-i . -i 
L j+i 




1 




r-i 


■ t\ ■ r 2 ■ h 




1 




r% 


■ h ■ r 2 ■ h 








r% 


■ h ■ r 2 ■ h 




h 






i ■ t\ 




k 


i 




i ■ t 2 




k 


i 




i ■ ts 




h 


i 




i ■ r 3 




T3 


-l 




i ■ r 2 




Tl 


■ i 




i ■ p 




P ■ 


i 



(CI) 



Here z — e l7T / 4 ,w = (1 — 2i)/y/E, diag(-) is the diagonal 
matrix with entries •, and permutations are specified in 
the standard way by the result of permuting the integers 
1 ... 8. Finally, X should be understood as acting after 
complex conjugation of the monopole field vector 



T:(C,0-A(2)P(Z)(C*,O- 



APPENDIX E: k-TENSOR 



(D7) 



The k [in Eq. I|4()[l] tensor was hiding the following 
expression 



APPENDIX D: GROUND STATE MANIFOLD 
PERMUTATIVE REPRESENTATION OF PSG 



6)4 = J2 K ij kl (kti 

ijkl 

= CiCiCo*Co* + 
+*C3&Ci*£o*- 
+<Ci&Ci*&*- 
-CiCoCo*Ci*- 

+«Co6C2*Cl* - 

+C2^oCo*6* -* 
-C2C1 Ci*C 2 * + 

+C3 ^3 C2* C2* — 

-C3C0O/C3* -i 
+C36 Ci*6* + 
-C36 C2* C3* - i 



C*C* 



C2 C2 Co* Co* + 
" C2 C3 Cl* Co ~ 
C 2 ClC3*Co* + 

6 6 Co* ft* - 
-CsCsCi*^*- 

C2 Co C3 Cl * + 

C3 Ci Co C2* + 
CoC 3 Ci*C 2 * + 

iClCoC3*C2* - 

C2 Ci Co C3 + 
Co C2 Ci* C3* + 

Co Co C3* C3* - 



*C3C3Co*Co 

~ C3 Ci C2* Co 

*Cl C2C3*Co 
«C2C3Co*Cl 

- C3 Co C2 Ci 

«CoC2 C3*Cl 

Ci C3 Co* C2* 
Co Co C2* C2* 

- Co Ci C3* C2 
CiC 2 Co*C 3 * 
Ci Co C2* C3* 
Ci CiC 3 *C 3 * 



*-CoCi Ci*Co 

* + Co C2 C2 Co 

*-CoC 3 C 3 *Co 

* + CoCoCi*Ci 
*-Ci C 2 C 2 *Ci 

* + Ci C3 C3 Ci 

- «C3CoCl*C2 
+ ClC2*C2 
*-C2C 3 C3*C2 

- «C2CoCl*C3 
+ «CoClC2*C3 

+ C2 C2 C3* C3* 
(El) 



In the basis of 1)38(1 the PSG realizes a permutative 
representation^ i.e. each generator can be written as 
G = A(G)P(G), where A(G) is a diagonal matrix with 
unimodular complex entries, and P(G) is a permutation 
matrix, acting on the 8-component vector (C,C)- O nc 
finds: 



(Dl) 
(D2) 
(D3) 
(D4) 
(D5) 
(D6) 



A(T0 


= diag(i,i, 1,-1, 1,1, 


P(Ti) 


= (43218765), 


A(T 2 ) 


= diag(z, z, iz, iz, — iz, — iz, z, z), 




= (21436587), 


A(^ (2) ) 


= diag(z ,z,—z ,—z,z,z ,z ,z), 




= (21437856), 


A(£(3)) 


= diag(2, 1, — z*, — 1, z, 1, z*, 1), 


^(£(3)) 


= (32417685), 


A(V) 


= diag(l, 1,1, -1,1, 1,1,-1), 


P(V) 


= (56781234), 


A(J) 


= diag(u;, w, iw, — w, — w, w, iw, w), 


P(J) 


= (86754231). 



APPENDIX F: PROOF THAT THE LOWEST 
DIMENSION OF A REPRESENTATION OF THE 
PSG IS 8 

We start by assuming that we are in a basis where T\ 
is diagonal. Denote one eigenvector by 



Tiipo = Xtpo- 



(Fl) 



The PSG algebra rule (|36|l dictates that T% cannot be 
diagonal in this basis. Using 

f 2 -f 1 =f 1 -f 2 (-i) , (F2) 

we find ipo is connected to a 4-cycle of eigenvectors 

Vm=T 2 m Vo, (F3) 

with eigenvalues 

f^ m = \i m i> m . (F4) 

This proves that the translations must all be constructed 
only of 4-cycles. As a consequence, any representation 
of the PSG can only be of a dimension d = 4 x n. 
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We know already that an 8 dimensional representation 
(rep) exists EQ - lB . 

so we need only consider the d = 4 
case, and show that it is impossible. 

Assume we have a d = 4 rep. Therefore, each eigen- 
value is non-degenerate, with a unique eigenvector. From 
(|Clj) we know that the inversion and T\ commute. It is 
easy to show that 

(i^ m )Az m = fi(iVm) (F5) 

follows, and since only a unique eigenvector has this 
eigenvalue, we conclude that 

(±Vm) = i>m (F6) 

and that in this rep 1=1. From the PSG algebra we 



now have 

in m = a-^i => n (3) = . (F7) 

The rotation 7t(3) is 3-fold, and therefore in this rep it 
is also unity! Finally, as 7Z(3) = 1 using the rule 

T 2 = £ (3 )Ti7^ (F8) 

we find T\ = T2, and we have a contradiction. 

This contradiction proves a 4-dim rep of the PSG is 
not possible, leaving us with the d = 8 as the lowest 
dimension of a rep of the PSG. 



1 L. Balents, L. Bartosch, A. Burkov, S. Sachdev, and 
K. Sengupta, Phys. Rev. B 71, 144508 (2005). 

2 L. Balents, L. Bartosch, A. Burkov, S. Sachdev, and 
K. Sengupta, Phys. Rev. B 71, 144509 (2005). 

3 C. Lannert, M. P. A.Fisher, and T. Senthil, Phys. Rev. B 

63, 134510 (2001). 

4 S. Sachdev and K. Park, Ann. Phys. (N.Y.) 298, 58 (2002). 

5 R. Moessner and J. T. Chalker, Phys. Rev. Lett. 80, 2929 
(1998). 

6 M. Hermele, M. P. A. Fisher, and L. Balents, Phys. Rev. 
B 69, 064404 (2004). 

7 P. W. Anderson, Science 235, 1196 (1987). 

8 X. G. Wen, Phys. Rev. B 65, 165113 (2002). 

9 R. Moessner and S. L. Sondhi, Phys. Rev. B 68, 184512 
(2003). 

10 D. Foerster, H. B. Nielsen, and M. Ninomiya, Phys. Lett. 
B 94, 135 (1980). 

11 O. I. Motrunich and T. Senthil, Phys. Rev. B 71, 125102 
(2005). 

12 J.-S. Bernier, Y.-J. Kao, and Y. B. Kim, Unpublished, 
|cond-mat/0411658| 

13 D. L. Bergman, R. Shindou, G. A. Fiete, and L. Balents, 
Unpublished, cond-mat/0510202 

14 H. Ueda, H. A. Katori, H. Mitamura, T. Goto, and H. Tak- 
agi, Phys. Rev. Lett. 94, 047202 (2005). 

15 H. Ueda, H. Mitamura, T. Goto, and H. Takagi, unpub- 
lished. 

16 K. Penc, N. Shannon, and H. Shiba, Phys. Rev. Lett. 93, 
197203 (2004). 

17 G. Misguich, D. Serban, and V. Pasquier, Phys. Rev. Lett. 
89, 137202 (2002). 

18 G. Misguich, V. Pasquier, F. Mila, and C. Lhuillier, Phys. 
Rev. B 71, 184424 (2005). 

19 A. Ralko, M. Ferrero, F. Becca, D. Ivanov, and F. Mila, 
Unpublished, cond-mat/0502294 

20 R. Moessner, S. L. Sondhi, and P. Chandra, Phys. Rev. B 

64, 144416 (2001). 



21 R. Moessner and S. L. Sondhi, Phys. Rev. Lett. 86, 1881 
(2001). 

22 D. A. Huse, W. Krauth, R. Moessner, and S. L. Sondhi, 
Phys. Rev. Lett. 91, 167004 (2003). 

23 D. S. Rokhsar and S. A. Kivelson, Phys. Rev. Lett. 61, 
2376 (1988). 

24 C. L. Henley, Phys. Rev. Lett. 62, 2056 (1989). 

25 C. L. Henley, Can. J. Phys. 79, 1207 (2001). 

26 J. Villain, R. Bidaux, J. Carton, and R. Conte, J. Physique 
41, 1263 (1980). 

27 J. Kogut, Rev. Mod. Phys. 51, 659 (1979). 

28 S. V. Isakov, K. Gregor, R. Moessner, and S. L. Sondhi, 
Phys. Rev. Lett. 93, 167204 (2004). 

29 L. Balents, L. Bartosch, A. Burkov, S. Sachdev, and 
K. Sengupta, Competing orders and non-landau- ginzburg- 
wilson criticality in (hose) mott transitions, pedagogi- 
cal article written for proceedings "Physics of Strongly 
Correlated Electron Systems" YKIS2004 workshop, held 
at Yukawa Institute, Kyoto, Japan, November 2004, 
|cond-mat /0504692 

30 E. Kroumova, M. I. Aroyo, J. M. Perez Mato, A. Kirov, C. 
Capillas, S. Ivantchev and H. Wondratschek. "Bilbao Crys- 
tallographic Server: useful databases and tools for phase 
transitions studies". Phase Transitions (2003), 76, Nos. 1- 
2, 155-170. 

31 B. I. Halperin, T. C. Lubensky, and S. K. Ma, Phys. Rev. 
Lett. 32, 292 (1974). 

32 E. Brezin, J. C. L. Guillou, and J. Zinn- Justin, Phys. Rev. 
B 10, 89 (1974). 

33 C. Dasgupta and B. I. Halperin, Phys. Rev. Lett. 47, 1556 
(1981). 

34 O. I. Motrunich and A. Vishwanath, Phys. Rev. B 70, 
075104 (2004). 



D. L. Bergman et al., unpublished, to appear. 



36 

36 The PSG is useful as a means of classifying emergent gauge 
structures as discussed by X. G. Wen in Ref. |j 



